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1 . GENERAL DESCRIPTION OF THE COMPUTER PROGRAM 


This computer program is used to implement the Cate-Liebig algorithm for con- 
tinuous data. The Cate-Liebig approach is intended to provide an alternative 
to stepwise multiple linear regression for the development of predictive 
models. In the program, cases are allocated to individual variables in an 
iterative manner based on the zero means of optimum binary (0-1) dummy vari- 
ables created by using the Cate-Nelson critical level algorithm. The means 
of the allocated cases, along with a calculated overall maximum applied at the 
final critical level, are utilized to define a linear response surface with an 
independent equation for each variable. 

The predicted value for any case is the minimum of all the values predicted 
for that case by the various individual equations. Variables are deleted from 
the model if less than a user-specified number of cases is allocated to that 
variable. Provision is made for testing the model on an independent set of 
data. 

This documentation and the program structure are modeled after the widely known 
biomedical (BMD) computer programs developed by the University of California. 

It is intended to be readily usable for persons familiar with the BMD program 
package; however, this documentation may be used without additional reference 
to any other publication. 



2. USER'S GUIDE 


2.1 PROGRAM MODES 

The program may be run in two primary modes. In the first mode, the program 
fits a set of dependent data using the Cate-Liebig algorithm and optionally 
provides a test of the resulting model on an independent data set. In the 
second mode, the program runs a user-specified model on an independent data 
set. 

2.2 PROGRAM OUTPUT 

Output from this program includes the following: 

1. Summary of allocations. 

2. Model coefficients. 

3. For fit and test data (if any): 

a. Limiting variable for each case. 

b. Estimate, actual value, and error of dependent variable. 

c. Simple error statistics, including mean error, mean square error, 
root mean square error, and coefficient of determination. 

d. Plot of estimated versus actual value of the dependent variable. 

4. Optional output, such as 

a. Plot of raw data for each variable. 

b. Summary of consistency check. 

c. Summary of critical-level splits after each iteration. 

d. Plot of cases allocated to each variable. 

2.3 LIMITATIONS PER PROBLEM 
Limitations per problem are as follows: 

1. Number of original variables (NOV); 1 < NOV < 80 



2. Total number of variables after transgeneration (IP); 1 _< IP £ 80 

3. Number of variable format cards (NVFC); 1 £ NVFC £ 10 

4. Number of cases in fit (N); 2 £ N £ 2000 

5. Number of test cases (NTC); 1 £ NTC £ 2000 

6. Number of transgeneration cards (NTGC); 0 £ NTGC £ 99 
This program allows transgeneration of variables. 

2.4 ORDER OF CARDS IN JOB DECK 

Cards in the job deck should appear in the following order (NOTE: items b 

through g may be repeated as often as desired): 

a. Program deck (a listing of the complete Fortran program is given in 
appendix A. All of the following cards are data read by that program.) 

b. Problem card 

c. Transgeneration cards (optional) 

d. "Labels" cards (optional) 

e. F-type variable format card(s) 

f. Data deck (place data input Jeck here, including test cases) 

g. Control -delete cards (optional). Alternatively when run in the model 
test mode only, the model definition cards are input at this point. 

h. Finish card 

2.5 CARD PREPARATION 


System cards (to be determined) 

Problem card: 


Column 

Problem 

1-6 

Problem (mandatory) 

10-15 

Alphanumeric problem name 

17-20 

Number of cases 
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Col umn 

Description 

24-25 

Number- of original variables 

29-30 

Number of transgeneration cards 

34-35 

Number of variables added by tra ns generation 

39-40 

Index of dependent variable 

44-45 

Number of control -delete cards 

49-50 

Number of labeled variables 

52-53 

Minimum R^ on split for allocation 

54-55 

Do not plot raw data if NPS = 1 

57-58 

Consistency check will not be used if = 1 

59-60 

Print summary of splits for each split if = 1 

64-65 

Do not print and plot each allocated variable if = 1 

66-70 

Number of test cases 

71-72 

Number of variable format cards 

73-75 

Smallest number of observations allowed for a variable 
to be allocated 

76-78 

Index of weighting vector 

NOTE: 

Any input variable may be used. All cases are assumed 
to be weighted equally if index is not given a value. 

If the number of cases = 0, the program is run in the model test mode and 

the model definition cards are required. 

Transgeneration card; 

Column 

Description 

1-6 

TRNGEN 

7-9 

Variable index K 

10-11 

Code from transgeneration list 

12-14 

Variable index I 

15-20 

Variable index J or constant C 


The term transgeneration is used to include transformations of input vari- 
ables and creation of new variables prior to computation. Notations to 
be used in the transgeneration list include the following. 
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, k = variable indices (need not be different) 
a constant 


a 2* a 3’*** = constants 

the number of cases, or sample size 

Transgeneration 

✓XT+ X k 

✓KT + /jq + 1 -> X ' k 

lo g 10 x i - X k 
e x l * X k 

arc sin ✓'XT’ -> X k 
arc sin JTTJTji + 1) 

+ arc sin v / 'nq~^n77TiT : rTJ •> X k 

v x i - x k 


x i + c - X k 


X i c * X k 
X i C ^ X k 


X i + x^ -> X 


x. - X, + X 

* vJ 


k 

k 


Vj 


x i /x j * x k 

If X, > c, 1 » X k ; 
otherwise 0 ->■ X, 


Restriction 

Xi > 0 


0 < X. < 1 
0 < (X./n) < 1 

x i t 0 

None 

None 

X. > 0 

None 

None 

None 

V° 

None 





Code 

Transgeneration 

Restriction 

16 

If X. > Xj, 1 -*■ X^; otherwise 0 -* X k 

None 

17 

l°9 e X| * X k 

X. > 0 

18 

If X. < c, 1 -* X k ; otherwise 0 X^ 

None 

20 

sin X. -»■ X^ 

None 

21 

cos X . **■ X^ 

None 

22 

arc tan X. -*• X k 

None 

23 

x/j * x k 

X. > 0 

24 

c X i » x k 

c > 0 

When a violation of a restriction in the right-hand column occurs durii 
transgeneration, the program will print a diagnostic message. 

Label 

cards: 


Col umn 

Description 


1-6 

Labels (mandatory) 


7-10 

Number of the variable (or category or 
(this number must be right-justified) 

index) to be named 

11-16 

Corresponding alphanumeric name 


17-20 

Number of another variable 


21-26 

Corresponding alphanumeric name * 


67-70 

Number of another variable 


71-76 

Corresponding alphanumeric name of that 

variable (up to 


7 variables per card) 
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Label cards allow the user to substitute alphanumeric names for the usual 
numeric indices that appear on the printed output. There may be from one 
to seven pairs of variable numbers and labels on each ''labels" card. If 
desired, only one pair may be specified on each card. However, the total 
number of labels appearing on all the "labels" cards must be equal the 
number of labels specified on the problem card. It is not necessary to 
label all the variables. Those labeled may be listed in any order. 

e. F-type variable format cards: The format is keypunched beginning with a 

left parenthesis, followed by a sequence of specifications, and closed 
by a right parenthesis. The first two specifications below are followed by a 
comma, except when they precede a slash or a right-hand parenthesis. 

Blank columns within the format are ignored, and columns 73 to 80 are not 
used. (A complete description of formats can be found in any Fortran 
programming manuals.) The specifications for F-type variable format cards 
are as follows: 

• "nFw.d" is the floating point (decimal ) indicator; n is the number of 
fields of width w (includes a sign and decimal point if punched); and 
d is the number of digits to the right of the decimal point if the 
decimal is not punched (0 _< d < w). If the decimal is punched, d is 
ignored. If n is not specified, it is assumed to be 1. 

• "mX” (alphabetic X), in which X is the skip indicator and m is the 
number of columns to be skipped. 

• "/" (slash) indicates, "go to the next data card." Depending on its 
location in the format statement, the "/" will either direct the pro- 
gram to go immediately to the next card (ignoring any further informa- 
tion on the current card) or skip one card altogether. For example, 
if a format begins with the program will automatically skip the 
first card, read the second, skip the third card, etc. If a format 
ends with "/," the program will automatically read the first card, 
skip the second, read the third card, etc. 

• "//" indicates, "go to the card after next." Two slashes "//" will 
direct the program to skip two cards. Any number of slashes may be 
used. 
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f. 


Control -delete cards 


Col umn Description 

1-6 C0NDEL 

7 Control value for first variable 

8 Control value for second variable 

• 

72 Control value for 66th variable 

If there are more than 66 variables, continue on another card of the same 
form until all variables have been specified. The variable numbers above 
refer to variables after transgeneration. 

NOTE: The following control values are permissible: 

0 - Linear form will be assumed. 

1 - Delete variable (or dependent variable). 

2 - A 0-1 form will be assumed. 

If no control -delete cards are included or if a field is left blank on 
the control -delete cards included in the deck, the value 0 will be 
assigned if the variable is not the dependent variable and the value 1 
will be assigned if it is the dependent variable. 

g. Model definition cards: If the number of cases on the problem card is 

set equal to zero, the control delete card(s) is (are) replaced by Card 1 

Column Description 

1-3 Number of variables in the model (integer) 

4-14 Maximum Y response in model (floating point) 

This card is followed by one card for each variable in the model in the 
following format. 

Col umn Description 

1-3 Variable number after transgeneration (integer) 
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Col unin Descri ption 

4-14 Intercept (linear model) or lower mean 0-1 model (floating 

point) 

15-25 Slope (floating point) 

26-36 Critical level (floating point) 

37-40 Type of response function for this variable (±1, 0) 

If the code for the type of response is set equal to zero, a linear model 
of the A + BX- will be assumed. If the slope, B, is positive. 


Y. + A + BX.j X.. <_ critical level 


Y.j = maximum otherwise 
If the slope is negative, 

Y i = A + BX.. X.. > critical level 


Y. = maximum otherwise 

If the code for the type of response is set equal to ±1 , a 0-1 mean model 
will be assumed. If the code equal +1, 

Y.. = lower mean X.. <_ critical level 


Y. = maximum 

otherwise 

If the code equals -1, 


Y. = lower mean 

X. > critical level 

Y. = maximum 

otherwise 


Note that the program output tabulates all of these values in the model sum- 
mary. When a 0-1 form has been assumed in the model fit, the appropriate 
lower mean value will be listed as the intercept. 


The form of the response functions may be mixed on a variable-by- variable 
basis. 


2-8 


i 




r- 

Er 


g. 


Finish card 
Col umn 


1-6 


Description 

Finish 


This card will notify the program that the entire job is finished. 
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APPENDIX A 


PROGRAM LISTING 

listing of the complete Fortran program is given on the following pages. 
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FILE . . : 


CL0516 FORTRAN PI 


C 

c 

c • • • 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

C. 

c 

c 


CAT 00060 
CAT00070 
CAT00080 
CATOOOOO 
CATO01O0 
CATOOUO 


NLV 

R2 

NPS 


CAT 00010 
CAT 00020 

• ..CAT00030 

THIS PHOGrtAM IMPLEMENTS THE CATE LIEBIG CRITICAL LEVEL ALGORITHM CATOOO'H) 

CATOOOSO 

PROBLEM CARD FORMAT 

COL NAME 

1-6 XMAN PRODLM 

10-1S CODE ALPHANUMERIC PROBLEM NAME 

N NUMMFR OF CASES 

NOV • NUMBER OF ORIGINAL VARIABLES 

NTGC NUMBER OF TRANSGENERATION CAROS CATOOUO 

NVA NUMBER OF VARIABLES AOOFD BY TPANSGENERATION CAT00130 

I NOFX I NOL X OF DEPENDENT VARIABLE’ CATOOUO 

NCO NUMHFR OF CONTROL DELETE CAROS CAT00150 

NUMBER OF LARELt'O VARIABLES CATOOUO 

VARIABLES WITH USD LESS THAN R2 WILL NOT CAT00170 

BE ALLOCATED CATOOUO 

00 NOT Pl.OI PAW DATA IF NPS = 1 CATOOUO 

NCON CONSISTENCY CHECK WILL NOT RE USED IF NCON = 1 CAT00200 

NS PRINT SUMMARY OF SPLITS FOR EACH SPLIT IF NS = I CAT002I0 

NPA 00 NOT PRINT AND PLOT EACH ALLOCATED VA&IAOLE IF NPa CAT00220 

NTC NUMBER OF TEST CASES CAT00240 

nvfc number of variable format cards catoo250 

NAL SMALLEST NUMBER OF OBSERVATIONS ALLOWED FOR A VARI ABLCAT00260 
. BE ALLOCATED CAT00270 

IDVTGINDEX OF WEIGHTING VECTOR CAT00200 

• •• ♦ ....CAT00290 

CAT00300 
CAT00310 
CAT 00320 
CAT00330 
CAT 00340 
CAT00350 
CAT 0Q360 
CAT 00370 
CAT00360 
CAT00390 
CAT00400 
CAT00410 
CAT 00420 
CAT00430 
CAT00440 
CAT 00450 
CAT00450 


1 7-20 
24-2S 
29-30 
34-35 
39-40 
44*45 
4 9-50 
52-53 

54-55 

57-58 

59-60 

64-65 

66-70 

71-72 

73-75 

76-76 


DIMENSION KTRANSO* 100) t TRANS <10.01 »PF.S (200 !»Y (0601 % XX (2000*3) 
1 • CU60) *ALBELU6Q) *AL!<U0*9) *CC U60 1 * AL ( 160* 10 ) 

DIMENSION XMI (?) ,XMA (?) *YPLOT (4000) *M (20001 *X (050 0.00) 

COM ON NT GC t TRANS * K TH ANS . Y t M , N I NCS t AL I * AL , C« XM I * XMA . YPLOT 
INTEGER C*ALDf L*CC *D 
DATA PR0BIM/4HOP0H/ *> 

DATA TPNCiF N/4HTRNG/ 

DATA CONDEl. /4HC0ND/ 

DATA FINI-SH/4HF INI/ 

NST = 0 


READ PROGRAM CARD 


C 
C 

c 

1 

100 " FORMAT ( A 4 i 5X » 2 X ♦ A 4 * 1 X * I *i * 6 { 3X * 1 2 ) t IX #F2.2* 12* 1 X 1 21 2 * 3X * 12* IS* 12* 21 CAT00470 


READ (5, 100) XMAN* CODE * N.NOV *NTGC *NVA* 1 NDEX *NCD * NLV *R? * NPS* NCON 
1 . MS , UP A t N r C f MVF C , f I AL . I D‘W TO 


l 3) 


2 

9001 

8020 

9999 

6021 

900? 

802? 

90 OA 
80?4 


8003 

7013 
7 Cl ) 

701? 

7005 

7008 


IF tXMAM-FlNISHI 2*9001*2 
I F ( AN - PRODLM 19002*3*9002 
Wf? X TE (6*6020 ) 

FORMAT (///?4H FINISH CARD ENCOUNTERED) 

WP I TE (6.8021 1 

FORMAT (l-RH PROGRAM TERMINATED) 

CALL EXIT 
WRITE (6. 002?) 

FORMAT ( 4 3 H NEITHER PROBLM NOR FINISH CARO ENCOUNTERED) 

GO TO 9999 
WRITF 16.8024) 

FORMAT 07H CARO INCORRECTLY PUNCHED OR MISSING.) 

GO 10 9999 
tPsNOV«HVft 

WRnf:(6*6003> CODE > N.NOV* InD£X*NVFC 

FORMAT (37H1 CATE LIEBIG CRITICAL LEVEL ANALYST S//2X* 1 2HPR0GPAM COOCAT0O64O 
IF* 17X* A6/2X* 1HHMUMCEH OF C A OF S * 1 6X * 1 4 / 2 X * 26HNUMBER OF ORIGINAL VARCAT 00050 
?IAMLf S.bX* 12 /2X .2-7H INDEX OF OFPl NDF.NT V AR I ABLE *<>X * I 2/2X *C AT 00660 

5 '3IHNUMUC4 Of VARIABLE FORMAT CAROS * 2X* 12// ) CAT00670 

t F (NT C . NF • 0 ) WR I T E ( 6 * 7 0 1 3 ) NTC CAT00600 

FORMAT <?X*21HMJ>*nFN OF TEST CASES *15) CAT00690 

IF(NOV.NE.O) W R I T F ( 6 • 7 0 1 1 ) NVA*IP CAT00700 

“ J*2X* CAT 0071 0 

CA T 00720 


CAT004G0 
C A T 0 0 4 9 0 
CATOOSOO 
CAT00510 
CAT 00520 
CAT 00530 
CAT 00540 
CAT005S4 
CAT00560 
CAT00570 
CAT 00560 
CAT 00590 
CAT 00600 
CAT 0061 0 
CAT 00620 
CAT 00630 


format r?x* AsMNuruuV of variables added by transgeneration , 12 * 
I2GMT0TAL NUMBER OF VARIABLES .13) 


IF (NLV.NF.O) WRITE (6# 70 12) 

FORMAT <?X»?1HVARIAHLKS ARE LABELED) 

IF ( NCR.NI. • 0) WRj 1 F { *» >7 005) 

FORMAT l ?X t 2SMCDMTROL DELETE CARDS USEO ) 
IF (NPS. HE. i > WRITE (6.7008) 

FORMAT (2X.13HPLOT PAW DATA ) 


CAT 00 730 
CAT00740 
CAT 00750 
CAT00760 
C AT 0 0 7 7 0 
CAT00780 


AinvnO UOM AO 
si aovd ivumnio 


r y*?* r XTK ■' - r - ■■ - 


**-3W*»» .. “WtKlT': 


filf.* 

7006 

7004 

7009 

9011 

901? 

7010 
9003 

80?3 

7 

7014 

C 

C . 

c 

600? 

9006 

6026 

7777 


9777 

6777 


977 * 

8775 

1000 

7015 

C 

c 

c 

106 

C 

C 

c 


c 

c 

c 


c 

c 

c 

12 

11 


WRITE { 6 • 90 1 1 ) 
LIMITED 


R2 


TO VARIABLES WITH A RSO GREATER 


. . . CL0516 FORTRAN R1 

IF(NPA.NE.l) WRITE (6*70081 

FORMAT (?X*30HPRINT AND PLOT EACH ALLOCATION ) 

IFCNS.FQ.l) WRITE (6 1 7 004) 

‘ FORMAT <2Xf41HPRINT SUMMARY OF SPLIT FOR EACH ITERATION) 

IF (NAL *NE . 0 ) WRJTF (6*7009) NAL 

FORMAT <2X,47HALL0CATI0NS LIMITED TO VARIABLES WITH AT LEAST 
1 t 1 2HORSERVAT IONS ) 

Ir (R2 • GE . 0 • 00 1 ) 

FORMAT (?X,57HALLOCAT IONS 
1HAN t FA • ?) 

IF (NOON.EQ* I ) WR I TE (6 • 901 2) 

FORMAT (2X*26MC0NSISTENCY CHECK NOT USED ) 

IF (IDWTG.NE.O) WRITE (6*7010) IOWTG 

FORMAT (?X*9HVAR1 ABLE • 1 4 , ?X , 26HUSED AS A WEIGHTING FACTOR ) 
IFUT.E0.1) WRITE (6*9003) 

FORMfiT <?X»?3MDATA INPUT IS FROM TAPE) 

IF (NOV «LE • 60 • AND. £P#LE*80) GO TO 7 
WRITE <6*«023) 

FORMAT (39H0MUMBER OF VARIABLES OUTSIDE OF LIMITS. ) 

GO TO 9999 

IF(NTGC«EO«0) GO TO 103 
WR1TE(6*7PIM 

FORMAT (//2Xt21HTRANSGENERATI0NS USED //2Xtl8H K 
00 1000 1-1 tttTGC 


CODE I J OR C ) 


READ TRANSGENERATION CARO 

READ (5*800?) YM AN. (KTRANS ( J, I ) , J=1 ,3) * TRANS <I> 

FORMAT <A4*2X*I3« 12* I3*F6.0> 

IF(YMAN-TRMGEN) 9006*7777,9006 
WRITE (6,6026) 

FORMAT (16H0TRANSGENERAT ION) 

GO TO 9004 
KIP = IP 

IF(NVA.LT.O) KIP = NOV 

IF(KT«ANS(1*I) .GT.KIP.0R.KTRANS(3,l).GT.KIP) GO TO 9777 
l F (KTRANS (2*1) .GT • 24 ) GO TO 9775 
GO TO 1000 
VRITE(6*«777> I 

FORMAT (///51H VARIABLE NUMBER SPECIFIED ON TRANSGF.NERAT ION 
14,1 ?H NOT CORRECT ) 

GO TO 9999 

WRITE (6*6775) KTRANS(?*I) 

FORMAT (///29H ILLEGAL TRANSGENERATION CODE ,I3,10H SPECIFIED) 
GO TO 9999 

WR JTE (6*7015) (KTRANS ( J* I ), J = l,3)fTRANS(I) 

FORMAT (3 (2X, 12) *2X,F7.2> 

READ LABELS CARO 

CALL R0LBL2 (NLV, IP, ALBEL) 

READ VARIABLE FORMAT CAROS 


IF(NVFC.LE.IO) GO TO 106 
WRI TF (6*4000) 

4000 FORMAT (1 H023X7 1 HNUMBER OF VARIABLE FORMAT CARDS INCORRECTLY 
X IED* ASSUMED TO BE 1.) 

NVFC - 1 
NVFC=NVFC*20 
READ (5, 6004) 

FORMAT ( 20A4 ) 

REWIND 2 
N3 = N ♦ NTC 
DO 112 K = 1 ,N3 
M (K) = 0 
NINCS=K 


106 

8004 


(RES ( I ) ,1 = 1 ,NVFC) 


READ DATA 

READ (5 , RES) ( Y (L ) *L s l v «NOV) 
1MNTGC) 11,11,12 

TRANSGFNERATION OF INPUT DATA 

CALL TRANGN 
CONTINUE 
WRITE (2,900) 'Y 


CAT00790 
CAT 00800 
CAT00810 
CAT 00820 
CAT00830 
• I 4, 2XC AT 00840 
CAT 0 0650 
CAT00060 
TCAT00870 
CAT 00880 
CAT00890 
CAT00900 
CAT 0091 0 
CAT00920 
CAT00930 
CAT 00940 
CAT00950 
CAT 0 0960 
CAT00970 
CAT00980 
CAT00990 
CATO 1000 
CATO 1010 
CAT01020 
CATO) 030 
CATO 1 040 
CATO 1050 
CATO 1 060 
CAT 01070 
CATO 1 060 
CAT 01 090 
CATO 1100 
CAT 0 1 1 1 0 
CAT 01120 
CAT 01 130 
CAT 01140 
CAT 0 1 150 
CATO 1 160 
CATO 1170 
CARO , I CAT 01180 
CAT C 1 1 90 
CAT 01200 
CATO 121 0 
CAT 0 1220 
CATO 1 230 
CATO 1240 
CATO 1 250 
CAT 0 1 260 
CATO 1 270 
CAT 01 280 
CAT 01 290 
CAT 0 1 300 
CAT01310 
CATO 1 320 
CAT 01 330 
CAT01340 
SPECIFCAT 0 1 350 
CAT 0 1 360 
CAT 01370 
CAT 0 1 380 
CAT 01390 
CAT 01400 
CA TO 1 4 1 0 
CATO 1420 
CATO 1430 
CA TO 1 440 
CATO 1450 
CAT 0 1 460 
CAT 0)470 
CAT 0 1 480 
CAT 01 490 
CAT 01500 
CAT 01510 
CAT 01520 
CAT 0 1 530 
CAT 0 1 540 
CAT 0)550 
CAT01S60 
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FILE, 

1129 

lit 1 


113 

360 


C 

C 

C 

35 If 
8011 

35? 

361 

9010 

8039 

36 
162 


37 

371 

C 

C 

C 


1221 


12? 


C 

C 

C 

1211 

i 21 

C 

c 


• . CL0516 FORTRAN PI 

00 1121 JJ a 1«IP 
X (K • JJ) = Y ( JJ) 

CONTINUE 
REWIND 2 
N1 a IP 

IF(N.EO.O) GO TO 213 

??<>W5 1 ’ 160 

cui=o 

NVAL = 0 

If (NCO.EQ.O) GO TO 36 
If (IP-66) 351 t 351 *352 

REAO CONTROL/DELETE CAROS 

RE AO (5 1 9011) WAN, (CCD .1 = 1* IP) 

FORMAT (A4.2X f 6611) 

GO TO 361 

REA0(5*80U) VMAN*(C(I)tI = lf66)*0UMt(C(I)tI=67tIP) 
IF ( VMAN-CUNOCL) 9010*36*9010 
WRITE (6*9030) 

FORMAT < 15MOCONTROL-OELETE) 

GO TO 9004 
C( INDEX) =1 

DO 162 I a It IP x 

C(I*IP) = 1 

TFCNAL.LT*!> NAL =1 

IFUOwTG.NE.O) C(IOWTG) a 1 

IQ = IP* 2 

00 371 I = 1*10 

00 371 J = 1*9 

ALI ( I « J) a. 0.0 

BEGIN SPLIT 

00 \A J=1#IP 

Iff C(J).EQ*1) GO TO 14 

0 a 0 

00 122 1*1 tN 
00 1221 K = 1*IP 
Y(K) a X ( I *K) 

IFfHCn.NE.O) GO TO 122 
0 = 0*1 
XX (0*3) a 1.0 

IFUOWTG.NE.O) XX (0* 3) = Y(IOWTG) 

YPLOT(D) = Y( INDEX) 

YPLOT (0*N) a Y ( J) 

XX (D*1)=Y(J) 

XX (0*2) =Y ( INDEX) 

CONTINUE 
N2 = 0 

IF(O.LT.NAL) GO TO 21 
!FCO*EQ. )) GO TO 21 
N3=N2 - 1 

IF CNPS.EQ. 1 ) GO TO 121 

PLOT RAW OATA 

WRITE <6t 1211 ) J* ALBELCJ) 

FORMAT C1M1*2?HRAW DATA FOR VARIABLE *I4*2X«A6/> 
CALL SETPLT(N) 

Kao 

SORT OATA INTO ASCENDING ORDER 

no 13 lal*N3 
A aXXUtl) 

A A=XX (1*2) 

AAA a XX ( I • 3) 

P =XX (1*1*1) 

BB=XX (1*1*2) 

PBB = XX ( I ♦ 1 • 3) 

IF (O.GE* A) GO TO 13 

XX(I.l) = 0 
XX ( I *2) aBB 
XX (1*3) a BBB 
XX (1*1*1) = A 


CAT01570 
CAT01560 
CATO 590 
CATO 1600 
CAT 01610 
CAT01620 
CATO 630 
CATO 1 640 
CAT 01 650 
CAT 01666 
CATO 167 0 
CAT016B0 
CAT01690 

m\m 

CAT01720 
CATO! 730 
CATO 1740 
CAT 01750 
CAT 0 760 
CATO 1770 
CAT 0 1700 
CATO l 790 
CAT 0 800 
CAT01810 
CATO 820 
CAT 0 1830 
CATO 640 
CAT 01650 
CATO I860 
CAT01B70 
CAT 01800 
CAT 0 690 
CAT 0 900 
CAT01910 
CAT 0 920 
CAT 0 930 
CATO 940 
C A T 0 1 95 0 
CATO i960 
CATO 1970 
CATO 960 
CATO 990 
CAT 02000 
CAT02010 
CAT 02020 
CAT02030 
CAT 02040 
CAT 02050 
CAT 02060 
CAT 02070 
CAT020H0 
CAT 02090 
CA TO? 1 00 
CAT 02110 
CAT 021 20 
CAT02130 
CAT 02140 
CAT02150 
CAT021 60 
CATO? 1 70 
CAT02160 
CAT02190 
CAT 02200 
CAT 0??1 0 
CAT 02220 
CAT 02?30 
CAT02240 
CAT 02250 
CAT 02260 
CAT 02270 
CAT02260 
CAT02290 
CAT02300 
CAT 0231 0 
CAT 02320 
CAT 02330 
CAT 02340 
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13 


CL 0516 FORTRAN PI 


135 


136 


137 


138 


7234 


7235 


XX (1*1.21 
XX(I*1.31 
CONTINUE 


=AA 


AAA 


GO TO 121 
j #N 3 


/! 


IF(K.EQ.l) 

88=1? 3 ,,a 
XA * 0.0 
XB * 0.0 
XBAR =0.0 
YBAR =0.0 
SUHXY =0.0 
SUMXSO =0.0 
SUMYSO =0.0 
Y1=0.0 
Y2=0.0 
Y1C = 0.0 
Y2C = 0.0 
NNN= NN*1 
XN=0 

lFCXXCNNtlJ.GE.XXCNNN.il) GO TO 133 
00 135 1= l.NN 
XN=XN*XX ( 1.3) 

XA = XA ♦ XX<I«n*XXUt3) 

Y1=Y1*XX(I.?1*XX(I»3! 

Y1C = Y 1C ♦ XX ( I . 3) 

YHAR=YRAK>XX ( I .21 «XX ( I .31 
00 136 1 =NNN.N2 
XN=XN*XX(If3! 

XBAR=X8AR*XX (1.3) 

XB = XB ♦ X X ( 1 . ] ) *XX (I. 3v 
Y2 = Y?*XX (It?) *XX (I .3) 

Y2C = Y2C ♦ XX ( I . 31 
Y9AR=YBAR*XX < I . 2) *XX (1.3! 

XBAR=XBAW/XN 

ybar=ybar/xn 

Yl=Yl/YiC .»• 

Y2=Y2/Y2C 
XA = XA/Y1C 
XB = XB/Y2C 
00 137 1= l.NN 

f 0. 0-XBAR1 '* (XX (1*2) -YBAR1 *XX (I «3! 
SUMXSQ=SUMXSO ♦ ( < 0. O-X0AR1 **21 *XX (1.3) 
SUMYSO=SUMYS(J ♦ ( (XX C I .2) -YBAR) **2) *XX (1.31 
00 1 35 I=NNN. N2 

SUMXY=SUHXY^* (1 .0-XBARl*(XX(It21-YBAR)*XX(I,3! 
5UHXS0=StlMXS0 ♦( () .0-XBAR1 **2> *XX ( I .31 
SUMYSO= SUMYSO + ( (XX ( I . 2) -YBAR) **2) *XX ( 1.31 
R = SORT (SUMXSO*SUMYSO) 

IF(R.OT.O.OOl) R = SUMXY/R 
RS0=R**2 

IF (RSO.LT. ALI ( J .71) GO TO 133 

AIKJ .7) =RSO 
ALI (J .21 =XX ( I I t 1 1 
AL I ( J .11=1.0 
ALT ( J. 9) = (XXdl + l.ll 
ALI ( J.fl) = XA 
1 1=3 
12=5 

IF(Y2«GT.Y1) GO TO 1382 
ALI (J . 1 1 =~1 .0 
ALKJ.9) = (XX(II*1.11 
“ = XB 


* XX ( I I . 1 ) 1 / 2.0 


133 

14 


103 

107 


t M1=Y1 
tI21=Y2 
•11*1 1 =NNN~1 
* 12*11 =N2-NN 


Arf , ( 6| xxm. ni/2.0 

ALI ( J.81 
11=5 
12=3 

1382 ALI ( J 
AL I ( J 
ALI (J 
ALI (J 
CONTINUE 
CONTINUE 

IF (NS.NE • 1 1 GO TO 168 
WRTTE (6. 1 031 

FORMAT (1 HI .4X.17HSUMMARY OF SPLITS//! 

WR I TF (6 1 1071 

FORMAT (//23X.flHCRlTICAL t5X .6H0 MEAN.5X .6HNWBER.9X • 6HNUM8ER/9X 
163MVARIABLE SIGN LEVEL X Y CASES 1 MEAN CASES 

? 00 15 1=1. IP 


CAT02350 
CAT 02360 
CAT 02370 
CAT02380 
CAT02390 
CAT02400 
CAT02410 
CAT02420 
CAT02430 
CAT02440 
CAT 02450 
CAT02460 
CAT 02470 
CAT02480 
CAT02490 
CAT02500 
CAT 02510 
CAT02520 
CAT02530 
CAT 02540 
CAT 02550 
CAT 02560 
CAT 02570 
CAT 02580' 
CAT02590 
CAT026C0 
CAT 02610 
CAT02620 
CAT02630 
CAT02640 
CAT02650 
CAT02660 
CAT02670 

CAT02700 
CAT027I0 
CAT02720 
CAT 02730 
CAT 02740 
CAT 02750 
CAT 02760 
CAT 02770 
CAT02780 
CAT 02790 
CAT02800 
CAT 0201 0 
CAT02820 
CAT02830 
CAT 02840 
CAT02850 
CAT 02660 
CAT 02070 
CAT02880 
CAT 02090 
CAT 02900 
CAT 0291 0 
CAT 02920 
CAT02930 
CAT 02940 
CAT02950 
CAT 02960 
CAT 02970 
CAT 02980 
CAT02990 
CAT03000 
CAT03010 
CAT 03020 
CAT 03030 
CAT 03040 
CAT 03050 
CAT03060 
CAT 03070 
CAT03080 
i CAT03090 
RSQCAT03100 
CAT03I 10 
CAT 03120 
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CL0516 FORTRAN PI 


\V 

c 


IF « C(I).EO.l) GO TO 15 A , 

•WRITE <6* 10a) ALRFLU) .I.(ALl(I*J) *J=1*2) t ALK 1*8) • < AL ^ 
FORMAT < 2X » A6. 3X. 1 3.F8,« 0*F8.2tF9.2*F8.2.F5.0#F9*2«f 6. 0* 
CONTINUE 


BEGIN ALLOCATION 


16* 


XM1N=99998.0 

IVARs-l 


sa s,e 


00 17 1= 1*10 

if < cin.r.o.n go to 17 


169 


I F C AL I (It7> • L T . R2 ) GO TO 17 
IF(CCU) .Nb'.O) GO TO 170 
IKCALl (!•«>. LT.NAL1 GO TO 17 


170 


GO TO 17 J 


171 


176 


IVA = ALIU.1) 

IFUVA.fO.CCU) ) GO TO 169 
IFCAII U.4).LT.NAl) GO TO 17 
IF(AL1(1*3) .GT.XMIN) CO TO 17 
I VAR = I ♦ IP 
GO TO 176 

IF (At I -C!»3> .GT.XMIN) GO TO IT 
1VAR=I 


17 


IV = I 

XMIN=AII (1*3) 

CON! I hue 

IFUVAR.EO.-l) GO TO 21 


C C T VAR ) = CUV> 


1791 


C - 


CCUVAR)=ALIUV*l) 

KVAR= I VAR 

IF ( I VAR .GT *Nl ) KVAR = IVAR - N1 
00 191 1 = 1 tN 
DO 1791 K = If IP 
Y LK ) - XUtK) 

IF (MU) .Nr. 0) GO TO 191 
IF(ALI(KVAR,1) .LT.O.O) GO TO 183 


VARIABLE ALLOCATED POSITIVE 


179 


IF I Y (KVAR) .GT.ALI (KVAR.2) ) GO TO 191 
MU) = IVAR 
GO TO 191 


C VARIABLE ALLOCATED NEGATIVE 
C 

183 


191 


20 


C 

C 

C 

21 

109 


|F(Y<XVAR) 

GO TO 179 
CONTINUE 
NVAL=NV4tU 
AL CNVAL? i 1 « IVAR 
DO 20 1 = 1.9 

AL (NVAL 1 1 ♦ 1) =AL 1 (KVAR • I ) 
NPS = 1 
GO TO 37 


.LE.ALI(KVARt2)) GO TO 191 


INITIAL ALLOCATION COMPLETE 


(I t J> »J=3t7)CAT03M0 
F5.?»I0F2.0)CAT03l50 
CAT03I60 


CAT03170 
CAT 03180 
CAT03190 


;AT03200 
CAT0321 0 
CAT03220 
CAT 03230 
CAT03240 
CAT 03250 
CAT 03260 
CAT 0327$ 
CAT03280 
CAT03290 
CAT 03300 
CAT0331 0 
CAT03320 
CAT 03330 
CAT03340 
CAT 0^350 
CAT 03360 
CAT03370 
CAT03380 
CAT03390 
CAT03400 


CAT03410 
1)3420 


22 


1092 


221 


WR I TE (6 • 109) 

FORMAT ( (HI *4X*19H INITIAL ALLOCATIONS) 

WR ITE(6«107) 

DO 2? 1=1 » NVAL 
L = ALU*)) 

J = L 

IF(L.GT.Nl) L-t L • N| 

WRITE (6* 104) ALBEL(L) * Jt ( AL < I t JJ> t JJ=2* 3) * AL « I *9) • (AL I 
1) 

WRITE (6* 1092) NAL .R2 

FORMAT (//■/■*> X .25HVARI ARLES WITH LESS THAN .I3.26M CASES 
1LOCATED. /5X • 45HVARI ARLES WITH SPLITS HAVING A RSO LESS 
1 ?0H WERE NOT ALLOCATED.) 

DO 221 I = ! * IQ 
DO 221 J =1*9 
ALT ( I * J) = 0.0 
00 222 I = 1 » NVAL 
L = ALU* I ) 

ALKL.9) = ALU. 10) 

AL I (L.5) = AL (NVAL *6) 


CAT0341 
CAT 03430 
CAT03440 
CAT03450 
CAT03460 
CAT03470 
CAT03A80 
CAT03490 
CAT03S00 
CAT03510 
CAT03520 
CAT03S30 
CAT 0 3540 
CAT03550 
CAT 03560 
CAT03570 
CAT03580 
CAT03590 
CA1 03600 
CAT03610 
CAT03620 
CAT03630 
CAT03640 
CAT 03650 
CAT 03660 
CAT 0 3670 
CAT03680 
CAT03690 
CAT 03700 
CAT 03710 
CAT03720 
CAT03730 
CAT 0.1 74 0 
CAT 03750 
CAT03760 
CAT 03770 
I*JJ) ♦ JJ=4» 8CAT 03780 
CAT03790 
CAT 03800 
WERE NOT ALC AT 0381 0 
THAN *F4.?.CAT03R?0 
CAT03830 
CAT 03840 
CAT 03850 
CAT03860 
CAT03870 
CAT 03880 
C AT03890 
C AT03900 
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ALI (L « 1 ) 
AL1 (L.2) 
AL I <L*3) 
15) ) 

ALTCLtA) 
15) I 

ALML.4) 


1 1 1 2) 
tr.3) 

. I (L 1 3) < 


ALI(U4) = ALI(Lt4) 
CONSISTENCY CHECK 


( AL I (L* 3) *ALI (L*4) 
(ALI CL tfl) # ALI (Lt4) 
ALI(UO ♦ AL(ItS) 


♦ AL I It 4) *ALf ItSI l/(ALI(Lt4) 

♦ AL < I #9) # AL tl« 5) )/(AL I (Lt4) 


CAT03910 
. CAT03920 
ALCI*CAT03930 


CAT03940 

ALU#CAT03950 


00 223} I * 
L = AL I ( 1 . 1 ) 


IFfL.ro. 0) CU> = 1 
IF (NCON.EQ. 1 ) GO TO 2136 
00 2220 L s ltNVAL 
2 = AL(L.l) 

IF (C ( I ) .EO. 2) GO TO 2120 


ALI. 7) = ULI iS) - AL I ( I .1) ) / ( ALI ( I • 9) - ALKI.8)) 

AL 1(1*6) = Al. I ( I «3) - ALI (I *7) *ALI (1*6) 

GO TO 2220 

AL! (1*7) a 0.0 

ALI ( I « 6) r ALI ( It 3) 

CONTINUE 
VR I TK (6*2219) 

FORMAT () HI , 4 X ♦ 1 7HC0NS I STENCY CHECK//54X . 8H0RIGINAL 1 20X » 3HNEW/5X * 
1 86HCASE t OR IGINAL NEW . ACTUAL Y LIMITING X EST Y 

2 LIMITING X EST Y ) 

00 2222 I * 1 *N 
XM1 = AL(NVALt6) 

JJ = 0 * x 

00 2232 K =^1 * IP 
Y (K) = X U t K ) 

DO 2221 J s itTO 
: IF (C < J) .EG. 1 ) GO TO 2221 

IF(J.GT.Nl) K r J - Nl’ 

XM 3 ALI < J.6I ♦ ALI (J.7)*Y(K> 

IF (Al.I f J. 1 ) .EO.-l • AND. Y (K) .LE.ALI ( J*2) > XM = AL (NVAL «6) 

IE < fi t M^ Vi 7: FO xi*Sl 4D i Y i? , * 6T *ALIIJt2)> XM = AL(NVALf6) 

IF(XM.GT.XMI) GO TO 2221 

XM 1 = XM 

JJ S J 

CONTINUE 

IF (MCI) .EQ.JJ) GO TO 2222 
J = M ( I ) 

IF(J.EO.O) GO TO 2222 
K=J 

IF(J.GT.Nl) K=J-N1 
XM2=ALI (J.6) ♦ALI (J.7)*Y(K> 

KK - JJ 

IF(JJ.GT.Ml) KK = JJ - Nl 

WRITE (6.2210) I t J» JJtY (INDEX) *Y(J) tXM2*Y(KK) tXHl 
FORMAT (5X.I4.4X. 16 .3X. iGtSXt FI 0*5*2 <5Xt2F10. 5)) 

A I 79 = 41. 1 ( J.4) -1 .0 

IF ( A I 79.LE • • 1 ) GO TO 2318 

ALI ( J. 3 ) - (ALI CJ»3)*ALI (J.4) - Y ( INOEX ) ) / (At I ( Jt4> - 1,0) 

ALI ( J ♦ 8) = (ALI (J.8)*ALI (J.4) - Y IK) 1 ✓ ( ALI ( J* 4) -1.0) 

ALI (J.4) 3 ALI (J.4) - 1.0 
GO TO 2119 
\ ALT (J.3)=0. 

ALI (J.G) =0. 

ALI(J.4.)=0. 

C(J)=1 

) AL I ( JJ .3) = (ALI (JJ«3) # ALI (JJ.4) ♦ Y( INDEX))/ (ALI ( JJ.4) ♦ 1.0) 

I F 7 JJ.GT .Nl ) K 3 JJ - Nl 

ALT (JJ.O) c (AU (JJ.8)*ALI(JJ f 4) ♦ Y (K) )/ (ALI (JJ.4 ) ♦ 1.0) 


ALI (JJ.4) 3 ALI (JJ.4) ♦ 1.0 
MO) = JJ 
CONTINUE 
WHITE (6.2225) 

FORMAT ( 1 H 1 . 4 X * 1 7HF INAL ALLOCATIONS ) 
WRITE (6.107) 

00 2226 I = It 10 

IF.7l-I.GT .Nl ) J s I - Nl 
IFIC(I) .EO.l) GO TO 2226 


CAT 03060 
CAT03970 
CAT 03900 
CAT03990 
CAT 04000 
CAT 040 1 0 
CAT 04020 
CAT 04030 
CAT 0404 0 
CAT04050 
CAT 04060 
CAT04070 
CAT04080 
CAT04090 
CAT 04 1 00 
CAT 04)10 


CAT 04120 
CAT 04 130 
CAT04140 
CAT04150 
CAT 04 1 60 
CAT04170 
CIT04100 
CAT04190 


CAT 04200 
CAT 042 10 
CAT04220 


CAT04230 
CAT04240 
CAT 04250 
CAT04260 
CAT04270 
CAT 04280 
CAT 04290 
CAT 04300 
CAT04310 
CAT 04320 
CAT 04 330 
CAT04340 
CAT04350 
CAT 04360 
CAT04370 
CAT 04 380 
CAT 04390 
CATO 4 400 
CAT 04410 
CAT04420 
CAT 04430 


CAT 04 44 0 
CAT 04450 
CAT04460 
CAT 0447 0 
CAT 04480 
CAT 04490 
CAT 04500 
CAT0451 0 
CAT 04520 
CAT 04530 
CAT04540 
CAT04550 
CAT 04560 
CAT 04570 
CATO450O 
CAT 04590 
CAT 04600 
CAT0461 0 
CAT 04620 
CAT 04630 
CAT04640 
CAT04650 
CAT04660 
CAT04670 


WRITE (At 104) ALHELIJ) . 1 . ( AL I ( I . JJ) . JJ = 1 .2) . ALI (I .8) . < At I < I . JJ) .CAT04660 
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CAT04690 
CAT047Q0 
CAT04710 
CAT04720 
CAT04730 
CAT 0*740 
CAT 04 750 
CAT04760 

.I3.2X.A6.2X »F5.0 »2X .5HCASES/CAT04770 

CAT04780 


1JJ-3.S) 

2226 CONTINUE 

2126 IFCNPA.KQ.IJ GO TO 212 
00 211 J = It 10 
lFCCCJI.EO.I) 60 TO 211 

IFCJ.GT.NI) JZ = J - N1 
WRITE <6. I 05) J.ALHEL ( JZ) .ALI (J.4) 

105 FORMAT 1 1H) «4X .20HVAKI ABLE ALLOCATED 
1) 

WRITE (6. 1051 ) 

1051 FORMAT C//SX.4MCASE *8X« I MX. 12X» 1HY) 

LI ~=°ALI < J»4) 

00 210 I = l.N 
00 2101 K = l.IP 
2101 YCK) - X(I.K) 

IFCMCl) .NE.J) GO TO 210 
WHITE C6» 1 05?) I.YCJZ) .Y(INOEX) 

1052 FORMAT CSX. 14 * 2F 13. 4.F20.2) 

L = L ♦ 1 

VPIOTIL) = YC INDEX) 

YPl.OT CL»N) = YCJZ) 

210 CONTINUE 

WHITE (6. 1 05) J.ALBEL) JZ) .ALI (J.4) 

CAl.l. SETPLT (1.1 1 

211 CONTINUE 
60 TO 212 

213 READ) 5.2 1 4 ) NVAL. AL (NVAL.6) 

214 FORMAT Cl 3. FI 1.5) 

10 = 2*IM 
00 215 I = 1. IQ 

215 CCI) = 1 

00 2)7 I = 1 .NVAL ' .. 

REA0C5.216) J.ALI (J.fl) .ALI ( J.9) .ALI ( J.2) .ALKJ.l) 

216 FORMAT < 1 3. 3F 1 1 .5.F4.0 L 

217 CCJ) = 0 

212 WRITE (6. 1 ) 0) 

110 FORMAT 1 1M1 »4X* 1 3HM00EL SUMHARY/5X . 55HVARI ABLE NUM 
1 SLOPE CR III CAL LEVEL ) 

DO 224 I = l.IQ 
L = I 

IFM.6T.Nt) L=I-N1 

IF C CCI). EO. U 00 TO 224 

IF (Fi.EQ.Q) 00 10 2242 

AU(!!9) I: = ‘ CM.I?ll5) 2 - 4 ALI Cl. 3) l/CALI Cl .9) - ALMI.8M 
ALICI. 81 = ALICI.3) - ALICI.9)*AL 1(1*0) 

00 TO 2242 

224) AL 1(1.9) = 0.0 

ALT (1.0) = ALI (1.31 

2242 CONTINUE . 

ViRI TF (6. Ill) ALBEL (L).I.AL ICI.8). ALI 1 1*9). ALI (1*2) 

224 CONTINUE _ . 

111 FORMAT C6X.A6.2X. 13.2 C2X.F1 1 .5) .6X.F11.5) 

WRITE (6. 1)4) AL (NVAL.6) 

114 FORMAT C//5X .7HMAX IMUM.Fl I .5) 

IFCN.EO.OI GO TO 21 A 
WR IT F. C 6.1001) 

1001 FORMAT C1H1.4X.14HDEPEN0ENT DATA //) 

CALL IfSTCN. INDEX. NVAL. NI.N5) 

210 IFCNTC.EO.O) GO TO 1 
WRITE (6. 1090). 

1 090 FORMAT C lNl . 4X.9HTEST-0ATA//) 

CALL TFSTCNTC. INDEX. NVAL. N1.N5) 

C.O TO I 

900 FORMAT C5E15. 61 

SUBROUTINE TESTCNTC. INOFX .NVAL »N1 *N5) . .. . ... 

OIMF.NSION All (J60.9) .AL 1)60.10) OHANSOOO) .KTRANS (3. 1 00) t V C 080) 
DIMENSION XMl (2) .XMAC2) .TPLOT (4000) «C ( 160) 

INTFGFR C 

COMMON iviniC* TPANStKTRANS* YoNtNINCS* AL 1 • AL«C» XMI , XHA t YPLOT 
VPTTF (6*10 1 ) 

101 FORMAT (SX«60MCASE VARlAOLE LIMITING X E5T Y ACT Y 

l ERROR! 


INTERCEPT 


CAT 04790 
CAT 04800 
CAT 0401 0 
CATG4620 
CAT04830 
CAT04840 
CAT04850 
CAT04060 
CAT 04870 
CAT04880 
CAT04A90 
CAT04900 
CAT 04910 
CAT 049?0 
CAT 04930 
CAT 04940 
CAT 04950 
CAT 04960 
CAT04970 
CAT 04980 
CAT04990 
CAT05000 
CAT 00010 
CAT 05020 
CAT 05030 
CAT 0504 0 
CAT050S0 
CAT 05060 
CAT 05070 
CAT 05080 
CAT C5090 
CAT 05100 
CAT05110 
CAT 051 20 
CAT 05 1 30 
CAT05140 
CAT 05150 
CAT 05160 
CAT05170 
CAT 051 GO 
CAT 05 1 90 
CAT 05200 
CAT052 1 0 
CAT05 ??0 
CAT 05230 
CAT05?40 
CAT05250 
CAT 05R60 
CAT05270 
CAT05280 
CAT05290 
CAT 05300 
CAT05310 
CAT 05320 
CAT05330 
CAT05340 
CAT 05350 
CAT05360 
CAT 05370 
CAT 05380 
CAT 05390 
CAT 054 00 
CAT 054 1 0 
CAT 05420 
CAT 05430 
CAT 05440 
CAT05450 
CAT0S460 
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220 

221 


22 


VBAR a 0.0 
REWIND 2 

IFINS.EO.il 60 TO 221 

N6=N5- 1 
00 220 I«ltN6 
REAO(2«900)Y 
CONTINUE 
00 22 I = ItNTC 
PE AD <2 .900) Y 
YHAR = YBAR ♦ Y I INDEX) 
VBAP = YHAR/NTC 
TSS a 0.0 


XRAR = 0.0 

= 0.0 


230 


RMSE 

REWIND 2 
IF(NS.EQ.l) 

oo ?30 I * 

PEAO(?.900) 


. _ j«?t 9 oo) y 

IMNS.EU.l) N 
f)0 2*1 1=1, NIC 


60 TO 231 
1 »N6 


CAT05470 
CAT05480 
CAT05490 
CAT 05500 
CAT0SS10 
CAT05S20 
CATOSf.30 
CAT0SS40 
CAT 05550 
CAT05560 
CAT0S570 
CATOS'iSO 
CAT0SS90 
CAT0S600 
CAT05610 
CAT05620 
CAT05630 


NIC 


READ(?*900) V 
XMIN=AC (NVAL,6) 

JJ= INDEX 
10 = 2*N1 
00 24 J a 1,10 
1FK(J),E0.I) 60 TO 24 

IFIJ.6T.N1) K a J - N1 
XM * ALI(J,A)> ALI(J,9)*Y(K) 


CAT0S660 
CATCS670 
CAT 05600 
CAT0S690 
CAT05700 
CAT0S710 


IF ( AL I ( J» 1 ) ,E0, 1 ,0,AN0» Y (K).GT.Atl(J,?)) XM = At CNVAL»6> 
IF(ALI(J»n.EO.- 1 . 0 .AND.Y{K).LE.ALi(J, 2 I) XM * AL(NVAC* 6 > 


24 


IF (XM.GT.XHIN) 60 TO 24 
XMIN = XM 

CONTINUE 

TSS = TSS ♦ (VRAM - YUN0EX)>**2 
fRROPsY ( INOEX) -XMIN ' 

XBAR = XBAR ♦ ERROP 

RMSE « RMSE ♦ ERROR’ERROR 

Z=Y(JJ) 

YPtOTin = V( INOEX) 


VPLOT ( I «N> = XMIN 
IFIJJ.EO. INOEX) “ 


1091 

241 


*99999,0 

XMIN, 


Y( INDEX) » ERROR 


1092 


102 


ARITE (6, 1091 • I,JJ*Z . 

FOPMAT(6X,I3,SX,I3»3C3X,EI0.4) ,F10.4) 

CONT 1 NUE 
XHAR s XBAR/NTC 
ESS * TSS - RMSE 
RSO * ESS/TSS 
XMSE * RMSE/NTC 
RMSE ” SORT (XMSE) 

WRITE (6.1092) XBAR.RMSF .XMSE.R5Q 

FORMAT <//5X. 10MMEAN ERROR . I 4X . F6. 2/5X .22HR00T MEAN SQUARE 
«F6.2/5Xt!7HKEAN SQUARE ERROR. 7X.F6.2/ _ M ^ 

' * . 5X»3HRSQ. 21X.F6.2) 

WRtTF (6. 102) 

FORMAT ( lH.lt SX.39HACTUAL Y ON Y-AXIS# ESTIMATE ON X-AXIS. 


(INI 

CALL_SETP1.T(NTC) 


2 

20 

3 


RETURN 
900 FORMAT (5E15.6) 

F MO 

SUBROUTINE POL RL2(NLBVAR.NVAR. ARRAY) 

INTEGER ALABEL. ARRAY. ITEST.DUMYtBLANK 
EQUIVALENCE (L ABEL t ALABf L ) 

DIMENSION ARRAY (60) . IDUM ( 7) .OUMY (7) 

DATA ALAUEL/3MLAB/ 

DATA BLANK/4H / 

DO I I=].NVAP 
APPAY ( I ) =BL ANK 
IF (NLBVAR) 9.9.2 
N=0 

RFAO (5.3) I TEST . ( IOUM( J) .OUMY(J) .J=1.7) 

FORMAT ( A3.3X * 7 ( 14 .A6) ) 

I E,iJ. TEST r l - AB EU. 4.6.4 
WRITE (6.5) 

FORMAT (36H0L ABELS CARD NOT FOUND WHEN EXPECTED) 
CALI. EXIT 
00 6 J=1 #7 


CAT 

CAT05740 
CAT05750 
CAT0S760 
CAT 05770 
CAT 05760 
CAT 05790 
CAT 05600 
CAT0SD10 
CAT05020 
CAT05630 
CAT05B40 
CAT 05B50 
CAT 05060 
CAT05870 
CAT058BO 
CAT 05890 
CAT 05900 
CAT 05910 
CAT 05920 
CAT 0H930 
CAT 05940 
CAT05950 
CAT05960 
CAT05970 
CAT 05960 
ERR0R.2XCAT05990 
CAT06000 
CAT06010 
CAT06020 
I CAT 06030 

CAT 06040 
CAT 06050 
CAT 06060 
CAT 06070 
CAT O60B0 
CAT06090 
CAT 06 100 
CAT06110 
CAT06120 
CAT 061 30 
CAT06140 
CAT08150 
CAT06160 
CAT 061 70 
CAT06180 
CAT 06 90 
CAT 06200 


CAT 06210 
CAT06220 


CAT 06230 
CAT Q6240 


A-l 1 




- 
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KsTOUH(J) 


CAT 06250 


IF <K) 11. A. 10 


CAT06260 

10 

fF(K-MVAM) 7.7.11 


CAT 06270 

11 

WRTTF. (6. 1 21 K.OUHY ( Jl 


CAT 06260 

12 

FORMA Til 8H0L ABELS CARO INDEX. 17* 18H INCORRECT. LABEL *A6.9H 

IGNORECAT06290 


10.1 


CAT06300 


GO TO 13 


CAT 063 10 

7 

ARRAY (Kl =DUMY (<JI 


CAT06320 

13 

N=N«1 


CAT 06330 

IFIN-NIBVAR1 6.9.9 


CAT06340 

8 

CONTINUE 


CAT 06350 


GO TO 20 


CAT 06360 

9 

RETURN 


CAT 06370 


FNIY 


’ CAT 06380 


SUBROUTINE TRANGN 


CAT06390 


DIMENSION TRANS (100) .K TRANS (3.100) 

«ALI 1160*9) *AL (160* 10) CAT06A00 


D IMF NS I ON X(OfiO). C ( 160) *XMI (?) *XMA (2) t YPLOT (40001 CAT06410 


COMMON NTGCt TRANS tKTRANStX.KtNINCSf ALI*AL*C. XMI *XMA* YPLOT 

CAT 06420 


INTEGER C 


CAT 06430 


AfAT <XX)=ATAN (XX/SORT <1.0-XX ## 2I> 


CAT 06440 


NCASE = K 


CAT06450 


00 100 IM.NTGC 


CAT06460 


M»KTRANS(1*I) 


CAT 064 7 0 


N=KTF?ANS<3.I> 


CAT06480 


NT RANS=K TRANS (2*1) 


CAT06490 

91 

l.rVH-Rl> 91.91.99 
TF(N-Gl) 92.92.99 


CAT 06500 
CAT06510 

92 

IF ( (NTP ANS-25 ) *NTRANS) 50*99*99 


CAT 06520 

99 

WRITE ( 6.199)1 

HISPUNCHEO OR OUT OF 

CAT 06530 

199 

FORMAT (22H TRANSGENERAT ION CARO *I3*27H 

0RDERCAT06540 


1) 


CAT 06550 


GO TO 100 


CAT06560 

50 

GO TO (1.2«3*4*5»6*7*8*9* 10*11* 12*13* 14* 
1.245.NTRANS 

15.16. 17. IS. 19,20. 21 

*22*23CAT06570 


* 

CAT 06580 

1 

I F ( X ( N ) ) 198.107*108 


CAT 06590 

107 

X<M)=0.0 


CAT06600 


GO TO 100 


CAT 06610 

108 

X f M ) =r SORT ( X (N) ) 


CAT 06620 

GO TO 100 


CAT 06630 

2 

IF ( X ( N ) ) 198.111.112 


CAT 06640 

111 

X(M)=1 # 0 


CAT 06650 


GO TO )00 


CAT 06660 

112 

XfMl-a SORT (X (N) ) *SORT <X(N)*1.0> 


CAT 06670 

GO TO 100 

* 

CAT06660 

3 

ZF(X(N>) 198. 190.114 


CAT 06690 

114 

X(M) r ALOGCXIN) >*0. 4342944019 


CAT067G0 

GO TO 100 


CAT 067 10 

4 

X CM) “EXP (X (N) ) 


CAT 06720 


GO TO 100 


CAT 06730 

5 

IF ( X (N) ) 198.107.117 


CAT 06740 

117 

IF(XCN)-l.O) 118.119*119 


CAT 06750 

118 

E=SORT ( X (N) ) 


CAT06760 

XIM|=ASNF<E) 


CAT 06770 


GO TO 100 


CAT 067 80 

119 

X (M) = 3.1415926536/2.0 


CAT06790 


GO TO 10G 

• 

CAT 06800 

6 

FNrNCASE 


CAT 068 1 0 


E=X(N)/(FN^1.0) 


CAT06620 


B=F ♦ 1 • 0/ CFN* 1 .0) 


CAT06030 


IF IF ) 198.123.124 


CAT06840 

123 

IF IB) 198.107.127 


CAT06850 

127 

X(M)=ASNF(SORT 18)) 

.. 

CAT06860 

GO TO 100 


CAT06870 

124 

IFCIU 198.128.129 


CAT 06880 

128 

X CM) =ASNF ( SORT (E) ) 


CAT08890 

GO TO 100 


CAT06900 

129 

e=^ort c e> 


CAT 069 10 

B=SOMT IB) 


CAT 06920 


X(M)=ASNF(E)^ASNF(8) 


CAT 06930 


GO TO 100 


CAT06940 

7 

IF IX IN)) 131.198*131 


CAT06950 

131 

X(M)r 1.0/XCN) 


C4T06960 

GO TO 100 


CAT 06970 

8 

X ( M ) = X ( N ) ♦ TRANS ( I ) 


CAT 06980 

GO TO 100 


CAT06990 

9 

XfMTrXCN)* TRANS ( I ) 


CAT 07000 


GO TO 100 


CATO 70 1 0 

10 

I F C X C N ) ) 198.107.133 


CAT07020 


A-12 


«x* % 

^moiho 


V 


r vK v ^ r fr > / a 
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133 

11 

12 

13 

14 

134 

15 

16 

17 

163 

16 

19 

20 

21 

fh 

187 

23 
168 

24 
189 

197 

198 
201 

100 


X (Ml * X IN) ••TRANS ( I) 
GO TO 
NEVRa 
X ( M 1 = 

GO TO 
NFWRa 


n 

100. 
T KA 


fPANS( I) 

X(N)#X(NF.VB) 

100 

TRANS ( 1) 

X(MUX(N)-X(NEWB) 

GO TO 100 
NF.Vft* TRANS ( II 
X(M)=X(N)*X(NEWB) 

GO TO JOO * 

NFWR= TRANS ( II 

IF (X (NF #B) ) 134* 197* 134 

XIM) = X(N)/X(NEWH) 

GO TO 100 

IF(X(N) .GT.TRANS(I) ) GO TO 111 
GO TO 107 
XNEWB=TRANS( I) 

1F(XCN>-(XNFWR>) 107# 111 till 
IF (X (N) ) 1 98 • 198 « 163 
X( M)s AI.OG(X(NI) 

GO TO 100 

IF (X IN) «LE. TP ANS (I ) ) GO TO 111 

X (M) =6.11*EXP( (•! .762042621 E 05*5. 59760 791 5E03*X 
1CN1 **21/ (1.254 162E05t273.0*XCNn) 

GO TO 100 

( X (N) ) 


TO 
X(M) = 
GO TO 
X(M)s 
GO TO 


100 

SIN 

100 

COS 

100 


( X < N) 1 


320 

330 


IF (X (N) -1 #57079632) 186*186# 198 
IF(X(N> #1.57079632) 193*187*187 
*X (M) =AT AN (X (N) ) 

,60 TO 100 
NEW8= TRANS ( I) 

IF (X (N) ) 198.198*188 
X<M)=X(N)**X(NEWB) 

GO TO 100 

IF (TRANS ( I) ) 198.198*189 
X f M ) = TRANS ( 1 ) **X (N) 

GO TO 100 
N=NF.WB 

WRITE t 6.201 )N.NINCS*KTRANS (2* I ) *K 

FORMAT (23M THE VALUE OF VARIABLE .14. 9H IN CASE 
1 THE RESTRICTIONS FOR TRANSGF.NERAT ION OF TYPE .13. 
2AM CONTINUED LEAVING VARIABLE tU.llH UNCHANGEO.) 
CONTINUE 
RETURN 
FNO 

SUBROUTINE SETPLT(N) 

DIMENSION KTRANS (3.100) * TRANS (100) * W < 080) . AL It 160 

01MI.NSI0N X**I <2) .XMA (2) .YPLOT (4000) .AL (160*10) t 

COMMON NTGC*TRANS*KTRANS.W*M*NINCS*ALI*AL.C* XMI 

INTEGER C 

NPV = N 

NC * 1 

K=1 

DO 330 1=1.2 
XMI ( I) =99999999.0 
XMA ( I) =-99999999.0 
DO 320 J=K.NPV 

XMIU)r AMIN1 (YPLOT (J). XMI (II) 

XMA ( I) = AMAX 1 (YPLOT (J) t XMA ( I ) > 

K=K#M 
NPV=NPV#M 
YMAX = XMA < 1 1 
YMIN=XMI(1) 

XMAX=XMA (2) 

XMIN=XMI(2) 

DO S50 1=1 *N 
2 = YPLOT (M#I) 


520 

530 

550 

570 

580 


Y = YPLOT (I) 

PLO 

INUC 

CALL PL0TR(Y*YHIN*YMAX9Z# 


CALL PL0TR<Y.YMIN*YMAX#2* 
CONTI 


xmin.xmax.no 

XMINfXMAX.NO 


CAT07030 
CAT 07040 

hwsh 

CAT 07070 
CAT07080 
CAT07090 
CAT 07100 
CAT 07110 
CAT07120 
CAT07130 
CAT 07140 
CAT07150 
CAT 07160 
CAT07170 
CAT 07180 
CAT07190 
CAT 07200 
CAT 0721 0 
CAT07229 
CAT 07230 
CAT 07240 
CAT 07250 

V CAT 07260 

(N)-2.050772636*XCAT07270 

. CAT07280 

CAT07290 
CAT07300 
CAT07310 
CAT 07320 
CAT07330 
CAT07340 
CAT07350 
CAT07360 
CAT07370 
CAT 07380 
CAT07390 
CAT07400 
CAT 074 1 0 
CAT07420 
CAT07430 
CAT 0744 0 
CAT07450 
CAT 07460 

* I5.55H VIOLATED CAT07470 

1H./40H THE PROGRC A T O' 7 4 8 0 
CAT07490 
CAT07500 
CAT07S10 
CAT 07520 
CAT 07530 

*9) *C ( 160) CAT 07540 

CAT 07550 

.XMA .YPLOT CAT 07560 

CAT 07570 
CAT07580 
CAT07590 
CAT07600 
CAT07610 
CAT07620 
CAT07630 
CAT07640 
CAT 07650 
CAT07660 
CAT07670 
CAT 07680 
CAT07690 
CAT 07700 
CAT077 1 0 
CAT 07 720 
CAT 07730 
CAT 07740 
CAT07750 
CAT07760 
CAT 07770 
CAT 07780 
CAT07790 
CAT 07800 
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8000 
8010 
8020 
80 30 
8040 


10 


20 

30 

40 

50 

60 


70 

80 

no 

100 

no 

120 


130 


140 

150 

160 

170 


180 

190 


200 

210 

220 


230 

240 

2S0 

260 

270 

280 


RETURN 
‘ END 

S IJMROUT VNE PLOTR (X » ?MIN*ZMAX * Y. WMIN* WMAX *NC) 

01 MENS TON Y (1 5) .CLABU2) *GF<10) •FMTCI21 tXY <51 « 101 ) * 5VM 115) 
INTEGER XY* BLANKS 
DATA SYM<1>/1H,/ 

NP e 1 

DATA TC*TPt8LANKS/lH.*lH*,lH / 

DATA OF/ 4H 1 X v 1 4 H 2X**4H 3X f 94M 4X»94H 5X**4H 

1 4H 7X • 9 4H 8X * ♦ 4H 9X..4H lOX/ 

DATA FMT/ • (17X»ft • * 1 5 (FI • * # 2 .3* • * » 8 X ) / • * »7X* •* ■ •iU(FP 
i • 8 X 1: : » • * * F 1 2 • 1 9 • 3 1 •'/ 

FOPPAT ( 1 H 6X5(F12.3*BX) *Fl?.3/l7X »5<F12. 3*0X1 ) 

FORMAT f 1 H F12*3*1X*103A1*F12*3) 

FORMAT OH 13X 9 103A1 1 
FORMAT 1 1H 1 4X 9 1 0 1 A 1 ) 

FORMAT U5X* 20 (5H^*..*) tlHO 
OATA NCC/?/ 

IF(NCC) 109200*10 
KL = 0 

CALL SCALE ( WMIN* WMAX * 1 00* 0* JY * YMIN* YMAX 9 Y I J) 

yr=ymax-ymin 

U=JY 

IF (J*(J-lon70»30*30 
IF (KL ) SO* 59-4-4-Q' 

WRITE (0*80401 
IF(KL) 170*170*50 
CI.AB(l) = YMIN 
DO 60 1 = 2*11 
CL AB(I)=CLAR(I-1) ♦YIJ 

WRITE (6*6000) (CLAB(I) 9 1 = 1*11*2) 9 (CL ABU) tJ=2* 10*2) 

IF (KL) 40*40*430 

IF( J-5) 90* 100*80 

J = J-5 

JYT=5-J 

CONTINUE 

FMT(?)=GF(JY) 

IF (KL) 120*120*140 
FMT ( 7) =GF ( JY) 

TT = JY 

TT = TT *YI J/l 0.0 
rLA 8 (l)= YMIN*TT 

00 130 1=2.10 

CLAH ( I > =CLAM( 1-1 ) ♦YlJ 

WRITE ( 6 * FMT') (CLAB (I) *1=2*10*2)* (CL AB(I) *1 = 1*9 *2) 

IF (KL ) 140* 140*430 

1 F ( JY-5 ) 150*160*150 

WRITE (6# 8030) (TC*I=1*J ) * (TP* (TC* 1 = 1 *4) *K=1 9 19) *TP* <TC* 

IF (KL) 170*170*120 
WRITE (6*8040) 

IF (KL) 170*170*120 
CONTINUE 
NCC = 0 
IC = 0 

IF (NP) 280* 180*180 
no 190 1=1*51 
DO 190 J=1*J01 
XY ( I* J) =BLANK5 

CALL SCALE (2MIN*ZMAX*50.0* JX*XMIN*XMAX*XIU) 

xr=xmax-xmin 

GO TO 200 
ENTRY PLOTS 
IF (NC) 4 20 * 360 *21 0 
IF (NP) ?H0 * 220*220 

00 270 N=1*NC 
SYMB=SYM(N) 

XOIFFR=XMAX-X 

IF (XOIFFR) 230*240*240 
XOIFFP = 0# 0 
VOIFFP=YHAX-Y(N) 

IF ( YD IFFP) 250*260*260 
YDIFFR = 0» 0 

1 =51 .0-(50.0*XDIFFR)/XR»,5 

K = 101.0-M0n.0*YDIFFR)/YR*.5 
CALL FORM?(SYMB*XYCL*K) ) 

CONTINUE 
GO TO 440 
00 290 1=1*101 


CATO701 0 
CAT07820 
CAT07830 
CAT07840 
CAT07850 
CAT07060 
CAT07870 
CAT07800 
6X* • CATO709O 
CAT07900 
••2 f 3***CAT07910 
CAT07920 
CAT07930 
CAT0794 0 
CAT 07950 
CAT07960 
CAT07970 
CAT 07980 
CAT07990 
CAT 08000 
CAT 08010 
CAT 08020 
CAT 08030 
CAT 0804 0 
CAT 00050 
CAT 08060 
CAT 08070 
CAT 08080 
CAT08090 
CAT00100 
CAT 08 1 10 
CAT 06120 
CAT 08 1 30 
CAT 08140 
CAT 00 1 50 
CAT 08 1 60 
CAT 0 8170 
C A T 0 0 1 6 0 
CAT08190 
CAT 00200 
CAT08210 
CAT 00220 
CAT 08230 
CAT 08240 
CAT 00 250 
CAT 00260 
CAT 08270 
I=1*JYT)CAT08200 
CAT 08290 
CAT 00300 
CAT 003 1 0 
CAT 00320 
CAT08330 
CAT 00340 
CAT0O350 
CAT00360 
CAT 06370 
CAT 00380 
CAT 00 390 
CAT 084 00 
CAT 084 1 0 
CAT00420 
CAT 08430 
CAT 00440 
CAT 08450 
CAT08460 
CAT 00470 
CAT 08480 
CAT0B490 
CAT00509 
CAT 0851 0 
CAT 08520 
CAT 08530 
CAT 00540 
CAT 08550 
CAT 08560 
CAT0R570 
CAT08500 
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XY^jtllsGLANKS 


300 

310 

320 


330 


340 

350 


360 


370 

380 


390 


400 

410 

4?0 


430 

440 


00 320 N=1*NC 
SYMR=SYM<N) 

YOIFFRsYMAX-YCN) 

K=101.0-<10n.0»YDIFFR)/YR* # 5 
CALL FORM? <SYMB*XV (L*K) ) 

IF <M0D(IC*5/ 1340*330*340 
W=TP 

GO TO 350 
W=TC 

WRITE (6*0O10>X*W*<XYU*N> *N=1*101)*W*X 

IC=IC*1 

GO TO 440 

M=6-JX 

LL=50*M 

f=JX 

IF (5-JX) 370*370*380 
T -0 • 0 

RLAB=XMAX-(T«XIJ)/5.0 

W=TC 

K=52 

DO 410 L=M*LL' 

K-K*j 

I=H0D(L»5> 

IF (I-l>400. 390.400 

SR!F„‘^«9H ,R *-AB,W,(Xy(K,N).N=l.l01).W.RUB 
RL A 8 -RLAo^X N 
W = TC 

GO TO 410 

coNTiNue ,8020)w * <XY(K ^ N) ,N:sl>lon » w 

KL =1 

GO TO 20 •> 

NCC=l 

PETURN 

END 

SUBROUTINE SCALE CYMIN*YMAX *YINT* JY* TYMIN*TYMAX t YI J> 


SUBROUTINE •SCALE 1 CALCULATES THE SCALING FOR •PLOTR* 

/l* 9* 1 *5*2* 0*3* 0*4 *0*5*0* 7#5* 10»0«1S*0*20*0/ 


10 


DIMENSION CUO) 

DATA C 

DATA TEST / 0.76293945E-05/ 
yr=ymax-yhin 

TT=YR/Y INT 

IF m.LE.0.0) TT = UO 
J = ALOGIO (TT) ♦TEST 
E=10.0**J 
TT=TT/E 


20 

30 

40 


1=0 

1F(“ 


50 

60 


70 


80 

90 


100 


no 

120 


(TT-l #0 4 TEST) 20*30*30 
TT=TT»10.0 
E=E/10.0 
I = 1 ♦ 1 

IF (9-1 ) 40*50*50 
E=E*10.0 

IF(TT-C(I) >60*70*30 
Y I J=C ( II # E 
GO TO 90 
Y=YMIN/CU) 

J=Y 

T 

IF (O.OOOl-ABS(T-Y)) 80*60*60 
YlJsCll*l) # E 

X = t (YMAX 4 YH1N)/YIJ-YINT >/2*0**00001 

IF(K> 100*120*120 
Y=K 

IF(X-Y) 110*120*110 

K=K-i 

TYMlNsK 

tymin=yij*tymin 

TYMAX=TYHIN*YINT*YIJ 


CAT08590 
CAT08600 
CAT08610 
CAT08620 
CAT 00630 
CATOB640 
CATOA650 
CAT 08660 
CAT 08670 
CAT 08680 
CAT 08690 
C4T08700 
CAT 087! 0 
CAT08720 
CAT0R730 
CAT08740 
CAT08750 
CAT 06760 
CAT08770 
CAT08760 
CAT08790 
CAT 08800 
CAT08810 
CAT08820 
CAT 08830 
CAT08840 
CAT 08850 
CAT 08860 
CAT08870 
CAT00880 
CAT 08890 
CAT 06900 
CAT 0891 0 
CAT 08920 
CAT 08930 
CAT08940 
CAT 08950 
CAT08960 
CAT08970 
CAT 08980 
CAT 08990 
CAT 09000 
CAT 090 1 0 
CAT 09020 
CAT 09030 
CAT09040 
CAT 09050 
CAT09060 
CAT09070 
CAT09080 
CAT 09090 
CAT09100 
CAT 09110 
CAT09120 
CAT09130 
CATC9140 
CAT 091 50 
CAT 09160 
CAT 09170 
CAT09180 
CAT 09190 
CAT 09200 
CAT 09210 
CAT 09220 
CAT 09230 
CAT 09240 
CAT 09250 
CAT 09260 
CAT09270 
CAT 09280 
CAT 09290 
CAT 09300 
CAT 09310 
CAT 09320 
CAT09330 
CAT09340 
CAT09350 
CAT 09360 
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... ( YMAX-TYM4X-TEST ) 130*130*30 CAT09370 

130 Y1JJ=C(!)«E . . CAT 093B0 

XT=( (YMAX*YMIN)/YIJJ-YlNT)/2.0*. 00001 CAT09390 

KT=XT CAT09400 

IF (KT ) 1 AO* 150* 150 CAT09410 

140 YT=KT „ CAT 09420 

. IF (XT.NE.VT) KT=KT-1 CAT09430 

150 T YMINT=KT CAT09440 

TYMINT=YI JJ»TYMINT CAT09450 

TYHAXT=TYMINT*Y1NT»YIJJ CAT09460 

IF (YMAX-TYMAXT.GT.TEST) GO TO 160 CAT09470 

TYMINsTYMINT CAT09480 

TYMAX=TYHAXT CAT09490 

YIJ=YIJJ CAT09500 

K=KT CAT09510 

160 TT =VINT/1 0.0 CAT09520 

.JY = TT*. 000001 CAT09530 

YI J=YINT« (YlJ/10.0) CAT 09540 

J=TYMIN/ YIJ CAT 09550 

IF (K) 170* IB0*1S0 CAT09560 

170 J=J-1 CAT09570 

180 J=J*JY*JY-K CAT095B0 

JY=J CAT09590 

PETURN CAT09600 

ENO CAT09610 

SUBROUTINE FORMS (SYMB.XY) CAT09620 

INTEGER XY«SYMO. BLANK/* */*TEST ( 18) /»2 *»*3 ».*4 ***5 ***6 ** CAT09630 

1*7 • * * B * * • 9 * . * A * * • R ***C • * *0 • *•£ • * *F *.»G * * *H * . * I **•/ »/ CAT09640 
IF (XY, EO, BLANK ) GO TO 20 CAT09650 

00 10 1=1*17 CAT09660 

IF (XY.NE.TE.ST ( I ) ) GO TO 10 CAT09670 

PUT IN NFXT SYMBOL OF ARRAY FOR MULTIPLE POINTS CAT09600 

XY=TEST (1*1) . CAT09690 

GO TO 30 < CAT09700 

10 CONTINUE • CAT097 1 0 

IF (XY.FO.TEST (18) ) GO TO 30 CAT09720 

’ IF OTHER THAN CHARACTERS IN ARRAY TEST PUT IN CHARACTER 2. CAT09730 

XY = TEST ( 1 ) CAT 09740 

GO TO 30 CAT 09750 

IF BLANK, PUT IN SYM30L CAT09760 

20 XY =SYMS CATC9770 

30 RETURN CAT 09700 

ENO CAT09790 
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SAMPLE INPUT DATA 
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APPENDIX B 


SAMPLE INPUT DATA 


Sample input data are as follows: 


FUF. . 

. 

L AMFI.S.0001 

<FP. 


01 

S 

10 

Op 

4 

10 

01 

s 

Q 

04 

f> 

H 

PS 

s 

Q 

ns 

c* 

10 

0/ 

4 

M 

0^ 

7 

10 

04 

S 

4 

10 

h 

Q 

11 

7 

O 

IP 

4 

A 

u 

s 

7 

14 

7 

7 

15 

7 

7 

is 

b 

A 

17 

4 

7 

1 * 

7 

7 

14 

s 

A 

?o <» 
rO'IDFL 
riMsn 

10 


M20TEST 

H?(j 20 

rji.K-iM (> 0 0 ? 

* ?f 5 • n t F*. 0 ) 
i P.5 

4 P.O 

? P.0 

5 1.0 

5 P.5 

i * . 0 

S P • 0. 

4 4.5 

1 1.0 

4 3.0 

J 3.0 
P P.o 

4 P.5 

5 3.S 

? <i • ft 

1 1.0 

3 P.O 

1 l.o 

P P.O 

1 1.0 


DATA 

ns oo 

H0003 


PI 

00 05 

C0004 


01 05.1000 1 1 

00005 Y 


0! 


01 


B-3 


A§£ U 


>Z^ 




20000 
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APPENDIX C 
SAMPLE OUTPUT DATA 



ra 


CATE LIEBIG CRITICAL LEVEL ANALYSIS 

PROGRAM CODE H20 
NUMBER OF CASES 20 
NUMBER OF ORIGINAL VARIABLES S 
INOFX OF DEPENDENT VARIABLE 5 
NUMBER OF VARIABLE FORMAT CAROS I 


NUMBER OF VARIABLES AOOEO BY TRANSGENERATION 
VARIABLES ARE LABELED 
CONTROL DELETE CAROS USEO 
PLOT RAM DATA 

PRINT AND PLOT EACH ALLOCATION 
PRINT SUMMARr OF SPLIT FOR EACH ITERATION 
ALLOCATIONS LIMITEO TO VARIABLES WITH AT LEAST 2 OBSERVATIONS 
ALLOCATIONS LIMITEO TO VARIABLES WITH A RSO GREATER THAN 0.10 


TOTAL NUMBER OF VARIABLES 






C-5 




' ’•* ?■ * **5#** *J*! X ■ ? !£*£:■' i' »•*■*- - * v.t- •••r Dr. 


RAW OATA FOR VARIABLE 1 NUMB 


4*000 6*000 12*000 16*000 20*000 

2.000 6.000 10.000 14*000 16.000 



;■■ iW-^:W4*w . .. ** i.v •. 4 f..aute, ... xm^A <smw*«. « «.*».: 




RAW DATA FOR VARIABLE 2 


4.000 


4.600 


5.600 


3*400 4.400 5.200 6.000 


0 

1 

OY 


4.125 


3.750 


3.375 


3.000 


2.625 


2.250 


1.675 


1.500 


1.125 


0.750 


0.375 



4.125 


3.750 


3.375 


3.000 


2.625 


2.250 


1.675 


1.500 


1.125 


0.750 


0.375 





8*800 9.600 10.600 

9.200 10.000 


8.600 


9.200 


9.600 


10.000 


f •••• 
10 . 60 ! 




RAW DATA FOR VARIABLE 4 
0.500 


‘Moo’ **** **** ***2.5oo* **** **** ***3.S5o* **** **** ***{.So8* 

2.000 3.000 4.000 





SUMMARY OF SPLITS 



CRITICAL 

0 MEAN 

NUMBER 


NUMBER 

GN 

LEVEL 

X 

Y 

CASES 

1 MEAN 

CASES RSO 

*1. 

15.00 

18.00 

1.40 

5. 

2.50 

15. 0.37 

1 . 

5.00 

4.50 

1.95 

10. 

2.50 

10. 0.12 

1* 

8.00 

7.50 

2.10 

10. 

2.35 

10. 0.03 

i. 

1.00 

1.00 

1.00 

4. 

2.53 

16. 0.61 




SUMMARY Of SPLITS 




NUMB VARI { BLE SI $ N 

H t 

£ 3 

0 4 


IvjN 

i: 


14.00 

5.00 

8.00 
2.00 


0 MEAN 
X Y 


7.50 

2.00 


NUMBER 

o c i? s 

3-J 1: 

2. 00 4 . 


. "EJN Til? RSO 

|»88 8 I olio 

i*71 li: S.-32 


WWSN*^** 

qf ooop. 0- 


C-10 




SUMMARY OF SPLITS 


NUMB 

H 

C 

0 


CRITICAL 

VARIABLE SIGN Lf ‘ 


1 : 1 $ 

3*00 


0 MEAN 

*!:» i:18 

1:88 1:88 


NUMBER 

CASES 


NUMBER 

MEAN CASES. RS0 


l: l:il *1: liti 


19 
2.75 


a# o«oi 


SUMMARY OF SPLITS 


VARIABLE SIGN 


CRITICAL 
LEVEL 


0 MEAN 


NUMB 

1 

1. 

11.00 

7.80 

3.10 

5. 

3.50 

1. 

H 

2 

1. 

6.00 

6.00 

3.00 

3. 

3.33 

3. 

C 

3 

•I. 

7.00 

0.20 

3.10 

5. 

3.50 

1. 

0 

4 

1. 

3.00 

3.00 

3.00 

2. 

3.25 

4. 


NUMBER NUMBER 

MEAN CASES RSO 
0.40 
0*50 
0.40 
0.25 


012 


, T , . 




SUMMARY OF SPLITS 


VARIABLE SIGN C (!lvEL AL 


1 : 




10.00 

6.00 

9.00 

4.00 


0 MEAN 

7.33 

6.00 

8.00 

5.00 


NUMBER 

3.17 

“.00 


3 

3.17 

3.25 


2 * 

3 . 

2 . 


NUMBER 
MEAN CASES RSO 
1. 0.33 


3.50 

3.50 

3.50 

3.25 


2 . 1.00 

1. 0.33 

2 . 0.0 







SUMMARY OF SPLITS 


VARIABLE SIGN 
•1 • 


CRITICAL 0 MEAN NUMBER „ „ NUMBER M)i 
LEVEL X Y CASES 1 MEAN CASES RSQ 

8.00 14.00 3.50 I. 3.50 I. 0.0 

0.0 0.0 0.0 0 . 0.0 0 . 0.0 

7.00 10.00 3.50 1. 3.50 1. 0.0 

4.00 5.00 3.50 I. 3.50 I. 0.0 



C-14 




. V v.l . 







CONSISTENCY CHECK 


CAf| ORIGINAL NEW 


ACTUAL Y 
2*00000 


ORIGINAS 


LIM|TJNG*X 


10000 


■ST Y 

f.00000 


NEW 


L, "m.5 ebb. 




' •vmhu fea 


■ TT-r . -M M r M rl , || „ , H M MW | BU | 


C— 1 7 


final allocations 



H 

0 


VARIABLE SIGN 

2 1 . 


4 1 . 


CRITICAL 

LEVEL 


0 MEAN 
X 

4.78 


1.78 


NUMBER 
Y CASES 
2.39 9. 

1*78 9. 









VARIABLE ALLOCATEO 


2 


H 


9. CASES 


X 

Y 

5.0000 

2.5000 

4.0000 

2.0000 

6.0000 

3.0000 

5.0000 

2.5000 

4.0000 

2.0000 

6.0000 

3.0000 

4.0000 

2.0000 

5.0000 

2.5000 

4.0000 

2.0000 



e-19 


If?'* W0& + 




VARIABLE ALLOCATED 2 
4.000 


9. CASES 


4.400 4*600 5.300 

4.300 4.600 5.000 5.400 


3.000 


2. 900 


2*800 


2.700 


2.600 


2.500 


2.400 


2.300 


2.200 


2.200 


2.000 


4 

4.000 


4.400 


4*800 


•S.zao* 


5.600 


•600 



C-20 



VARIABLE ALLOCATEO 




2.0000 

3.0000 

1.0000 


3.0000 


2.0000 


1.0000 

1.0000 

2.0000 

1.0000 



• & 















VARIABLE ALLOCATED 
1.000 

4 

1*200 

0 0* CASES 

1*400 

1*600 

1.800 

2.000 

2.200 

2.400 

2*600 

2.600 

3.000 


3.250 












3.250 

3.000 










2 


3.000 

2*750 












2.750 

2*500 












2*500 

2.250 












2.250 

2.000 





3 







2*000 

1.750 












1*750 

1*500 












1*500 

1*250 









■* 



1*250 

. 1.000 

♦ ‘ 











1*000 

0*750 










* 


0*750 

1*000 

1*200 

1 *400 

1*600 

1.800 

2.000 

»•£*•• . *^ 
2*200 

2*400 

’••SUlo**' 

2*800 

3.000 




^ J'JI. 1! !-l -l l w ^ ^ 


DEPENDENT OATA 



CASE 

VARIABLE 

LIMITING 

X 

| 

2 

5*0000 


2 

2 

4*0000 


3 

4 

2*0000 


4 

2 

6.0000 


5 

2 

5.0000 


6 

4 

3*0000 


7 

2 

4.0000 


8 

4 

4*0000 


9 

4 

1.0000 


10 

2 

6*0000 


11 

4 

3.0000 


12 

2 

4.0000 


13 

2 

5.0000 


14 

4 

5.00D0 


15 

4 

2.0000 


\h 

4 

i.0000 


17 

2 

4.0000 


ie 

4 

1.0000 


19 

4 

2.0000 


20 

4 

1.0000 


Me AN 

ERROR 


0.00 

ROOT 

MEAN SQUARE 

ERROR 

O.Ofl 

MEAN 

SQUARE error 

0.01 

RSQ 



0.99 


EST Y 

ACT Y 

ERROR 

2.5323 

2.5000 

•0.0323 

1*8871 

2.0000 

0.1129 

2*0000 

2.0000 

0.0„ 


3.0000 

2.5000 


3.0000 

3.0000 

0.0_ 

1.8871 

2.0000 

0.1129 

3.5000 

3.5000 

0.0 

1.0000 

1.0000 

0.0_ 

3.1774 

3.0000 

-0.1774 

3.0000 

3.0000 

0.0 

1*8871 

2.0000 

0.1129 

2.5323 

2.5000 

-0.0323 

3.5000 

3.5000 

0.0 

2.0000 

2.0000 

0.0 

1.0000 

1.0000 

0.0 

1.0871 

2.0000 

0*1129 

1.0000 

1.0000 

0*0 

2.0000 

2.0000 

0*0 

1.0000 

1.0000 

0*0 


C-23 


i , -,-rt HA 




C-24 


ACTUAL V ON Y-AXlSt ESTIMATE ON X-AXIS. 

1 .iroc ■ . 1 .SCO 2.400 3.000 3.600 

0.900 1.500 2*100 2.700 3*300 



i.^sa 

0.900 1*500 


’i.*«5* ; mu 

2.100 2.700 


FINISH CABO ENCOUNTERED 
PROGRAM TERMINATED 




